Multi-scale Modeling Approach to Acoustic Emission during Plastic Deformation 
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We address the long standing problem of the origin of acoustic emission commonly observed 
during plastic deformation. We propose a frame-work to deal with the widely separated time scales 
of collective dislocation dynamics and elastic degrees of freedom to explain the nature of acoustic 
emission observed during the Portevin-Le Chatelier effect. The Ananthakrishna model is used as it 
explains most generic features of the phenomenon. Our results show that while acoustic emission 
bursts correlated with stress drops are well separated for the type C serrations, these bursts merge 
to form nearly continuous acoustic signals with overriding bursts for the propagating type A bands. 

PACS numbers: 83.50.-v, *43.40.Le, 62.20.fq, 05.45.-a 



Acoustic emission (AE) is observed in an unusually 
large number of situations. For example, it is observed 
during crack nucleation and propagation in fracture of 
solids micro- fracturing process 0, martensite trans- 
formation 0, 13] , peeling of an adhesive tape [F, ^ , collec- 
tive dislocation motion etc 0,i. Clearly, sources that 
lead to AE signals in such widely different situations are 
system specific even as the general mechanism attributed 
to AE is the abrupt release of the stored strain energy. 
The phenomenon is used as a non-destructive tool in un- 
derstanding the sources and mechanisms generating AE. 

Acoustic emission during plastic deformation refers 
to high frequency transient elastic waves generated by 
abrupt motion of dislocations. AE studies in plastically 
deforming metals and alloys have been reported for over 
four decades [1]. The changes in the AE signals during 
deformation differs from one type of experiment to an- 
other and also on the sample type. Some correlations 
has been established between AE signals and the nature 
of stress-strain (cr — e) curves (gj-fllj . Conventional yield 
phenomenon is accompanied by a peak in AE pattern just 
beyond the elastic regime that decays for larger strains. 
In contrast distinct AE patterns are observed in the case 
of unstable plastic deformations such as the Liiders band 
and the different types bands in the Portevin - Le Chate- 
lier (PLC) effect |10l - [l3| . Such differences in AE patterns 
in different experimental conditions (and samples) can 
be attributed to the way dislocations respond to external 
forces. Theoretical approaches to AE are based on Green 
function approach that use specific model sources such as 
an expanding loop which generate AE . Clearly, such 
approaches cannot be useful if one is interested in fol- 
lowing the changes in AE occurring during the course of 
deformation since AE signals (as also stress) are averages 
over dislocation activity in the entire sample. Despite the 
vast literature on the subject, we are not aware of any 
model that predicts the nature of acoustic emission dur- 
ing the entire course of deformation. 

The purpose of the paper is to propose a theoretical 
frame-work to describe both dislocation dynamics and 
elastic degrees of freedom simultaneously since it is the 
abrupt motion of dislocations that transmits the kinetic 
energy to the surrounding elastic medium triggering the 



AE signals. We address the problem in the context of the 
PLC effect where the signature of the AE signals are well 
correlated with the types of serrations and band types 
observed at different strain rates p^| - [l3| . Our results 
show that for type C serrations, the AE bursts which 
are correlated with stress drops, are well separated. As 
strain rate is increased, the AE bursts tend to merge to 
form nearly continuous acoustic signals with overriding 
bursts for the propagating type A bands. 

It is well known that AE signals in the case of the 
Liiders and the PLC bands arise from collective behavior 
of dislocations [l0l-[l3l [Tsj . In such cases, it is necessary to 
simultaneously describe the collective behavior of dislo- 
cations and the elastic degrees of freedom. This so far has 
not been possible due to several difficulties. First, a ma- 
jor source of difficulty common to all plastic deformation 
experiments, is the absence of theoretical frame work to 
simultaneously treat the widely separated inertial time 
scale and that of dislocation dynamics. Second, there 
is lack of dislocation based models to describe collective 
behavior of dislocations [Tsj. Third, there is no clarity on 
how to describe transient acoustic waves. Finally, even 
in models describing collective dislocation motion, stress 
equilibration is assumed. This, however, no longer holds 
during the process of AE generation [l5l - [l7f . 

The PLC instability is characterized by three types of 
bands and the associated serrations [l^. On increasing 
strain rate or decreasing temperature, randomly nucle- 
ated static type C bands are seen, identified with large 
stress drops. Then the type B 'hopping' bands are seen. 
Here, a new band is formed ahead of the previous one 
in a spatially correlated way giving the visual impres- 
sion of hopping propagation. The serrations are more 
irregular with smaller amplitude compared to the type C 
serrations. Finally, the continuously propagating type A 
bands associated with small stress drops are seen. 

Our basic idea is to obtain the local plastic strain 
rate from model equations that describe the entire spatio- 
temporal evolution of plastic deformation and use it as a 
source term in the wave equation for the elastic strain. 
Here we use the Ananthakrishna (AK) model for the PLC 
effect [iB-Ell as it reproduces the band types [l5l4l7| , and 
several other generic features such as the existence of the 
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instability within a window of strain rates, the negative 
strain rate behavior etc [H, [l^. The model also pre- 
dicts chaotic stress drops which has been subsequently 
verified ^g^, HH . The basic idea of the model is that all 
the qualitative features of the PLC effect emerge from 
nonlinear interaction of a few dislocation populations, 
assumed to represent the collective degrees of freedom of 
the system. The model consists of densities of mobile, 
immobile, and decorated (Cottrell) type dislocations de- 
noted by pm(a;,T), pi,„(x,T) and Pc(x,t) respectivel y, in 
the scaled form. The scaled evolution equations are : 
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where t is the scaled time variable. The term 6oPm Eq. 
([1]), refers to the formation of dipoles and other disloca- 
tion locks, PmPim refcrs to the annihilation of a mobile 
dislocation with an immobile one and the source term pim 
represents the athermal or thermal reactivation of the im- 
mobile dislocation, apm represents the immobilization of 
mobile dislocations due to aggregation of solute atoms. 
Once a mobile dislocation starts acquiring solute atoms 
we regard it as Cottrell- type of dislocation pc- As more 
and more solute atoms aggregate, they eventually stop, 
and are considered as immobile dislocations pim- This is 
the source term apc in Eq. (j)^^jPm in Eq. ([T} rep- 
resents the rate of multiplication of dislocations due to 
cross slip. This depends on the velocity of mobile disloca- 

1 /2 

tions taken to be Vm{4>) = 4>Tf}-' where 0e// = {(jy—hp^j^ ) 
is the scaled effective stress, m the velocity exponent, 
and h a work hardening parameter. Further, cross-slip 
allows dislocations to spread into neighboring spatial lo- 
cations and thus gives rise to diffusive coupling (last term 
in Eq. ([T])). These equations are coupled to Eq. (jU that 
represents the constant strain rate deformation experi- 
ment. In Eq. ([4]), ia is the scaled applied strain rate, 
£{p,x,t) — pm{x,T)(f>™j j(x,t) is the local plastic strain 
rate, d the scaled effective modulus of the machine and 
the sample, and I the dimensionless length of the sam- 
ple. Note that Eq.Q assumes stress equilibration. The 
scaled constants, a, c and &o refer, respectively, to the 
concentration of solute atoms slowing down the mobile 
dislocation, the diffusion rate of solute atoms to mobile 
dislocations and the thermal and athermal reactivation of 
immobile dislocations. The relevant parameter is the ap- 
plied strain rate £a with respect to which different types 
of serrations and the associated bands are observed. The 
instability range is found in the interval 30 < < 1000. 
Equations (P-H]) are discretized on a grid of N points 



and solved using a adaptive step size differential equation 
solver ("MATLAB" 'odelSs'). In experiments, bands 
cannot propagate into the sample due to large strains 
at the grips. This is mimicked by choosing the bound- 
ary conditions pimi^jT) and pim{N,T) to be two orders 
higher than the rest of the sample. In addition, we im- 
pose pm(l,T) = Pc{N,t) = 0. The initial values of 
the dislocation densities are chosen to be uniformly dis- 
tributed with a Gaussian spread along the sample. For 
the numerical work, we use a = 0.8, & = 5 x lO^'', c = 
0.08, d = 6 X 10"^ TO = 3.0, h^O,D^ 0.25, N ^ 100. 

The above equations (Eqs. ([T]-!!])) are adequate to ob- 
tain the plastic strain rate only. However, noting that the 
abrupt collective dislocation motion triggers the transient 
elastic waves, we need to describe both elastic degrees of 
freedom and dislocation dynamics. This also implies that 
instantaneous stress following such an event will display 
fluctuations that damp-off in course of time. Indeed, the 
abrupt slip process induces dissipative forces that tend to 
oppose the accelerated motion of the slip interface. This 
is a mechanism that ensures eventual approach to me- 
chanical equilibrium. Following Ref. |22| . we represent 
this dissipation in terms of the Rayleigh dissipation func- 
tion (RDF) given by TZae — 2 
TZae with acoustic energy dissipated by noting that this 
has the form of the energy associated with abrupt disloca- 
tion motion during plastic deformation, i.e., TZae oc e^(r) 
[23} . Thus TZae is taken to be the energy of the tran- 
sient elastic waves. We have shown that the choice of 
representing the acoustic energy dissipated in terms of 
Rayleigh dissipation function has been successful in pre- 
dicting the nature of AE signals in varied situations such 
as the martensite transformation '3, '24*1, fracture [23I] and 
peeling of an adhesive tape 6, 25, 26). Writing down the 
kinetic energy J{e'^{y,t)dy, where p is the density), 
the potential energy (^ J e'^{y,t)dy, where p is the elas- 
tic constant), dispersion of the elastic waves (-yi^^]^ 
with D a constant), and dissipation TZae, we get (using 
Lagrange's equations motion). 
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(5) 
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The second and third terms (on the right hand side) arise 
from the dispersion and dissipation terms respectively. 
In addition, we have included the plastic strain rate (last 
term) calculated from Eqs. (HKHEl), and Eq. (g]). This 
acts as a source term in the wave equation for the elastic 
degrees of freedom that is expected to generate transient 
elastic waves. Note that Eq. ([5]) is general and applicable 
to any plastic deformation situation as long as the plastic 
strain rate is supplied. Transforming this equation into 
scaled variables used in the AK model, we have 
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the scaled and unsealed are 
= 6Vot where b is the Burgers 
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vector, /3, 7, 9 and Vq are constants used in the unsealed 
AK model equations. See Ref. [l^ for details.) 

Finally, appropriate boundary conditions needs to im- 
posed on Eq. ([6|) that should be consistent with those 
on Eqs. (HJU]). This however is not straightforward. To 
do this, we first note that numerical solution requires dis- 
cretization of Eqs. (lUIll) and Eq. (|6|). Further, as one end 
of the sample is fixed and a traction is applied to the other 
end, the total imposed strain rate is shared by the ma- 
chine and the sample. This implies that the machine elas- 
tic element should be included at both ends, i.e., the dis- 
crete form of the wave equation should contain equations 
of motion for the end points of the sample and machine. 
Then, the stiffness of the machine enters naturally in the 
equations for the end points. Then, boundary conditions 
of Eq. ([T]|3]) are automatically satisfied by these equa- 
tions. The relevant boundary conditions for discretized 
form of Eqs. ([6]) are ei(r) = 0, andejvlT) = iaT for r > 
where the subscript 1 and N refer to the end sites. The 
initial conditions are: £^(0) = + ^, i = 2, .., N — 1 with 
the random number ^ is drawn from interval — ^ < ^ < ^ . 

However, the time scale of plastic strain rate (i.e., 
Eqs.([l]|31)) is typically ~ while that of Eq. ^ is much 
smaller. Indeed, the step size in an adaptive step size 
algorithm used for the solution of Eqs. ([T]-|31) are signif- 
icantly larger that the time step required for integrating 
Eqs. (111). Thus, we need to ensure that the time vari- 
able in Eq. ^ and Eqs. (PU]) are mapped correctly. 
Denoting the i*'* integration time step in the AK model 
by ATi, for the time interval between r^+i < r < r^, 
we need to ensure that mSr' = Ati where St' is the 
fixed step size used for Eq. ([6]). Further, we use in- 
terpolated values for the plastic strain rate £p(fc,T) (for 
any fc*'' spatial element) obtained by using linear inter- 
polation formula Sp{k, t) — Sp{k, Ti) + '^p(*^''^^)~^p^^'"^'+i) t-^ 

where Ti < t < Ti^i where ti is z*'' time step of inte- 
gration of Eqs. (HJH]). Moreover, the plastic strain rate 
calculated from Eqs. ([T]|4]) has a much coarser length 
scale compared to the fine length scale required for wave 
propagation. Noting that the spatial coupling in the AK 
model appears only in Eq. ([T]), it is easy to show that the 
strain rate ep{x,T) in the AK model must be scaled by 
a factor (assumed to be constant) when used in the 
wave equation, i.e., ep{y,T) — \^ep{x,T) where x and y 
refer respectively to spatial coordinates in the AK model 
and Eq. (The range of A is 10^ - 10^.) The resuhs 
presented are for N = 100, A^ ~ 10"^, km = 5fcs, j^r^ — 
1500,7/(6lVo)2 = 10 and D/{eVaf = 1. Note that the 
velocity of acoustic waves is of right order for 9Vq ^ 100. 

Equations and Eq. ([6]) are in principle coupled 

since ip{y,T) is a function of stress 0(t). A self consis- 
tent solution of these equations is equivalent to solving 
the full dynamical problem involving both plastic defor- 
mation and elastic degrees of freedom with the attendant 
difficulties. This will not be attempted here. Instead we 
provide an approximate method akin to adiabatic meth- 





FIG. 1: (a) Uncorrelated type C bands for Sa ~ 40. (b) (Color 
online) Plots of stress and acoustic emission energy signals. 
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FIG. 2: (a) Partially propagating type B bands for ia = 130. 
(b) (Color online) Plots of stress and acoustic emission energy 
signals in asymptotic regime. The region between the arrows 
in figures (a) and (b) are identified. 
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FIG. 3: (a) Fully propagating type A bands at ia = 240. (b) 
(Color online) Plots of stress and acoustic emission energy 
signals in asymptotic regime. The region between the arrows 
in figures (a) and (b) are identified. 



ods. The procedure adopted is to first calculate ip{k,T) 
for the entire duration of time by solving Eqs. (IT]|4|). 
Then, Eqs. ^ is solved using ep{k,T) as a source term 
(along with the scale factor A). This gives the elastic 
strain ee{y,T). Then, the integral of ee{y,T) over the 
specimen dimension gives the transient stress (j)tr{T) ex- 
plicitly. While 4)tr (t) will be equal to (j) within the elastic 
limit, it will be different from (/) beyond this limit. 

The AK model predicts the three band types found 
with increasing strain rate [TB - flTj . At low ia, say ia — 
40, the uncorrelated static C bands are seen as shown in 
Fig. [TJa). The serrations are large and nearly regular. 
The scaled acoustic energy dissipated is obtained using 



dy 



dy, where F"* = T/{eVo)^. Since, 

stress drops in this case are due to isolated band nucle- 
ation, the AE pattern consists of well separated bursts 
that are well correlated with the stress drops. Figure 
[IJb) shows a typical stress-strain curve along with the 
AE bursts for ia = 40. The post burst AE is continuous 
that gradually increases until a new burst is seen jllj . 
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At intermediate strain rates, say ia = 130 hopping 
type B bands are seen as shown in Fig. [D^a). These prop- 
agate partially and stop mid-way. Another hopping band 
reappears in the neighborhood. Often, nucleation occurs 
at more than one location. The corresponding asymp- 
totic stress-time plot is shown in Fig. [2lb). The associ- 
ated serrations are irregular but are smaller in magnitude 
compared to the type C. While the correlation between 
stress drops and AE peaks still holds when the propaga- 
tion is short, the AE bursts are not as well separated as 
in the case of type C serrations. A plot of the AE signal 
is shown in Fig. Hfb). During hoping propagation, low 
level AE activity is seen in the region between two AE 
bursts (shown by arrows) (A few small bursts 

are also seen.) As we increase Ea, the extent of propaga- 
tion increases with concomitant decrease in stress drop 
magnitudes. At high we find fully propagating type A 
bands. Figure[3l^a) shows dislocation bands nucleating at 
one end of the sample and propagating continuously to 
other end for ia = 240. The corresponding AE pattern 
[Fig. [3] (b)] appears nearly continuous with a few over- 
riding bursts. Large bursts in AE are correlated with the 
nucleation of the band (or the band reaching the edge or 
due to occasional intersection of two bands). There is a 
low level AE activity during propagation (the region be- 
tween the arrows). In experiments, bands once nucleated 
trigger a burst in AE but during propagation very low ac- 
tivity is seen HI, [ll]- Thus, the generic features of 
AE signals during the PLC effect are well captured. 

In summary, we have developed a theoretical frame- 
work for dealing with widely separated inertial time scale 



and that of collective dislocation modes to explain the na- 
ture of acoustic emission patterns observed in the PLC ef- 
fect. This has been done by computing the plastic strain 
rate from the AK model for the PLC effect and using it as 
a source term in the wave equation. An important input 
in the theory is that the energy of the transient acous- 
tic wave dissipated caused by the abrupt slip (resulting 
from collective unpinning of dislocations) is represented 
in terms of the Rayleigh dissipation function 0, Q . The 
results show that for type C bands, well separated burst 
type AE signals that are correlated with stress drops are 
seen. As we increase the strain rate successive bursts tend 
to merge. For high Sa where type A propagating bands 
are seen, the bursts merge to form continuous type of AE 
signal. Over riding this are AE bursts that correspond 
to band nucleation or a band reaching the edge. These 
features are consistent with experimental results (TTl - [l3| . 
Other features such as those from hardening can not be 
captured here as there is very little hardening in the AK 
modelfH, [l^]- However, an extension of the AK model 
that removes this limitation can be used [131 ■ The frame- 
work is clearly applicable to other deformation conditions 
as long as a dislocation based model can be developed 
that captures major features of the phenomena. Finally, 
better approximate schemes have been designed that also 
give similar results [28|]. 
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